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ABSTRACT 


In  a  study  of  the  planetary -scale  vorticity  budget  of 
the  northern  summer  monsoon,  Holton  and  Colton  (1972)  found 
that  a  strong  vorticity  sink  is  needed  to  balance  the  gener- 
ation by  horizontal  divergence  at  200  mb.   In  this  work  we 
use  numerical  experiments  to  examine  the  hypothesis  that 
the  fluctuation  of  the  divergence  forcing  as  observed  by 
several  studies  serves  as  this  sink.   A  3-level  numerical 
model  is  used  with  the  thermal  forcing  specified  by  the  hori- 
zontal divergence  distribution  observed  by  Krishnamurti 
(1971).   Two  experiments  are  carried  out,  one  with  heating 
steady  in  time  and  the  other  fluctuating  in  time  with  a 
period  of  10  days.   The  time-mean  fields  of  the  quasi-steady 
state  solution  in  both  experiments  show  a  significant  west- 
ward phase  shift  of  the  Tibetan  high,  thus  the  fluctuation 
of  heating  cannot  be  a  sufficient  sink.   On  the  other  hand, 
the  positions  of  the  simulated  tropical  upper  tropospheric 
troughs  agree  quite  well  with  those  observed.   This  sug- 
gests that  damping  due  to  cumulus  transport  is  more  likely 
to  account  for  the  vorticity  sink.   Transient  wave  activi- 
ties are  also  found  at  the  upper  level  in  both  experiments. 
They  appear  to  be  drawing  energy  from  the  planetary-scale 
flow  and  could  serve  as  a  vorticity  sink.   These  waves  are 
most  active  in  the  fluctuating  heating  case,  suggesting 
that  local  barotropic  instability  at  upper  levels  may  be 
more  important  than  that  implied  by  the  time-mean  wind. 
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I.   INTRODUCTION 

In  the  study  of  the  northern  summer  planetary-scale  mon- 
soon circulation,  Holton  and  Colton  (1972)  discovered  that 
wery    strong  damping  is  necessary  to  balance  the  budget  of 
time-mean  planetary -scale  vorticity  at  the  upper  levels. 
Specifically  they  found  that  without  this  strong  damping  to 
counterbalance  the  generation  of  vorticity  by  horizontal 
divergence  the  amplitude  of  the  long  waves  of  the  monsoon 
circulation  would  be  too  large  at  the  200  mb  level  and  the 
phase  would  be  shifted  westward  by  up  to  5000  km. 

This  same  phase  problem  was  also  apparent  in  a  numerical 
study  by  Abbot  (1971),  who  used  a  2  layer  quasi -geostrophi c 
model  with  prescribed  heating  to  study  the  monsoon. 

In  looking  for  a  source  of  this  strong  damping  mechanism, 
Holton  and  Colton  suggested  two  hypotheses  which  may  explain 
this  strong  damping.   The  first  is  the  vertical  mixing  of 
vorticity  and  momentum  by  strong  cumulus  updrafts  and  down- 
drafts.   The  second  hypothesis  is  related  to  the  observed 
time  fluctuation  of  the  monsoon  system.   If  the  large-scale 
divergence  contains  a  strong  fluctuation  then  a  correlation 
term  of  the  divergence  and  vorticity  fluctuations  will  ap- 
pear in  the  time-mean  vorticity  equation  and  may  reduce  the 
generation  of  the  time-mean  vorticity  by  the  time-mean 
divergence. 
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Murakami  (1974)  used  a  three-component,  quas i -spectral 
model  with  prescribed  heating  which  has  a  fluctuating  compo- 
nent.  He  found  that  it  can  reasonably  simulate  the  ampli- 
tude and  the  phase  of  the  long  waves  associated  with  the 
northern  summer  monsoon  circulation.   But  his  model  suffers 
from  the  lack  of  the  zonal  basic  flow  and  the  resultant  ver- 
tical profile  of  the  divergence  field  differs  considerably 
from  observations.   Therefore  it  is  not  clear  which  mecha- 
nism causes  the  seemingly  correct  phase  and  amplitude  of 
his  model. 

The  purpose  of  this  study  is  to  examine  the  second  hypo- 
thesis of  Holton  and  Colton  using  the  non-linear,  finite 
difference  model  of  Monaco  and  Williams  (1975).   This  same 
model  was  used  to  simulate  the  northern  summer  monsoon  cir- 
culation by  Bellis  (1975).   In  his  study,  Bellis  found  that 
the  prescribed  heating,  which  was  specified  by  the  200  mb 
mean  horizontal  divergence  field  observed  by  Krishnamurti 
and  Rodgers  (1970),  can  reasonably  excite  several  important 
features  of  the  time-mean  monsoon  circulation.   However,  if 
extremely  large  values  of  damping  are  not  included,  a  phase 
shift  problem  still  exists,  although  the  amplitude  of  the 
disturbance  seemed  to  agree  with  that  observed. 

In  this  study  we  will  carry  out  two  controlled  experi- 
ments.  One  will  have  the  heating  specified  to  be  independ- 
ent of  time  and  the  other  will  have  a  fluctuating  heating 
component  similar  to  that  observed  and  used  by  Murakami. 
Hopefully  by  comparing  these  experiments  we  will  be  able  to 
evaluate  the  effects  of  this  fluctuation  in  the  monsoon  heating 
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Another  purpose  of  this  study  is  to  simulate  the  baro- 
tropic  instability  of  the  tropical  easterly  jet  associated 
with  the  northern  summer  monsoon.   This  instability  arises 
as  a  result  of  the  strong  easterly  jet  in  the  upper  tropo- 
sphere and  has  been  studied  by  Colton  (1974)  and  Bel  1  is 
(1975).   It  represents  an  important  scale-interaction  mecha 
nism  for  the  monsoon  circulation,  according  to  the  observa- 
tional study  by  Krishnamurti  (1971). 
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II.   THE  NUMERICAL  MODEL 

A.   BASIC  MODEL 

The  numerical  model  used  was  formulated  by  Monaco  and 
Williams  (1975)  and  is  similar  to  that  described  by  Arakawa 
and  Mintz  (1974)  who  used  the  primitive  equations  in  spheri- 
cal coordinates  with  "sigma"  as  the  vertical  coordinate. 
In  this  study  the  domain  of  calculation  is  restricted  to 
the  band  between  18°S  and  46°N  for  the  sake  of  computer 
economy,  and  a  free  slip  condition  is  applied  at  the  lateral 
boundaries.   The  model  has  been  previously  used  by  Bellis 
(1975)  for  the  area  of  0°E  to  180°E  only. 


B.   THE  GRID  SYSTEM 

The  horizontal  distribution  of  variables  is  illustrated 
in  Figure  1  which  is  the  so-called  C-scheme  arrangement 
(Arakawa  and  Mintz,  1974).   The  grid  spacing  is  4  degrees  in 
both  latitude  and  longitude.   The  troposphere  is  divided 
into  seven  equally  spaced  sigma  layers  in  the  vertical  with 
three  main  reporting  levels,  as  shown  in  Figure  2.   The 
sigma  coordinate  is  defined  as 

P  -  P  + 


a  = 


(2.1) 


where  p  is  the  pressure  of  the  sigma  level,  and  p.  is 
the  constant  pressure  height  of  the  tropopause  set  at  100 
mb,  and   it   is  the  terrain  pressure  defined  as 


TT  =   Ps   -   Pt 
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(2.2) 


The  surface  pressure,  p  ,  was  set  at  1013  mb  initially. 
The  vertical  boundary  conditions  are  -jr  =  0  at  o  =  0,  (p 
=  pt)  and  a  =  1 ,  (p  =  pg) . 

C.   SPECIFICATIONS  OF  HEATING  AND  DAMPING 

To  evaluate  the  effect  of  fluctuating  heating  as  sug- 
gested by  Holton  and  Colton  (1972),  two  controlled  experi- 
ments are  carried  out.   In  both  experiments  the  heating  func 
tion,  H,  contains  two  components: 


H  "  Hs  +  HAS 


(2.3) 


where   H   ,  the  zonally  symmetric  heating,  is  given  by  a 
continuous  adjustment  to  a  specified  equilibrium  temperature, 
T  ,  using  the  following  formula, 


=   rT  -  T  (y,q) 

T  -I 


(2.4) 


Here   T   is  the  temperature  and   x   is  a  specified  adjust- 
ment time  set  to  two  days  in  this  study.   The  equilibirum 
temperature,  T  ,  is  prescribed  to  maintain  the  zonal  mean 
north-south  temperature  gradient. 

The  horizontal  distribution  of  heating  is  derived  from 
the  200  mb  time-mean  divergence  field  observed  by  Krishna- 
murti  (1970)  as  follows: 


Q(°K  day"1)  =  Divergence  (sec-1)  x  106 


(2.5) 


The  field  is  smoothed  on  the  boundaries  by  reducing  the 
strength  of  the  heating  by  50%  at  the  points  just  inside 
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the  boundaries  and  80%  at  the  boundary.  The  entire  field 
was  smoothed  by  a  five  point  space  smoother  in  order  to  re- 
move small  scale  features.  In  the  vertical,  the  full  value 
of  heating  derived  by  (2.5)  is  specified  at  level  2  and  a 
vertical  distribution  is  assumed  such  that  the  heating  at 
level  1  and  3  are  50%  and  25%,  respectively,  of  the  value 
at  level  2. 

The  zonally  asymmetric  component,  H.~.  takes  the  follow' 


ing  form, 


2tt 


HA$  =  Q  x  (1  .0  +  A(SIN  $f-   *)) 


(2.6) 


where  A  is  the  amplitude  of  the  heating  fluctuation  whose 
period  is  T.   In  the  first  of  the  two  experiments  A  =  0;  in 
the  second,  A  =  0.5.   A  10  to  15  day  period  was  observed  by 
Murakami  (1974)  and  Krishnamurti  and  Bhlame  (1976)  in  the 
northern  summer  monsoon.   Here  ten  days  is  chosen  for  the 
periodicity  T  for  the  sake  of  computation  economy. 

Linear  damping  terms  are  included  in  the  momentum  equa- 
tions to  represent  the  effect  of  scale- independent  friction. 
Horizontal  diffusion  is  also  introduced  to  include  the  effect 
of  scale-dependent  friction  for  it  damps  the  small  scale 
fields,  which  are  generated  by  the  finite  difference  scheme 
in  the  model . 


D.   INTEGRATION  SCHEME  AND  INITIAL  CONDITIONS 

The  time  integration  is  comprised  of  continuous  sections 
of  one  Matsuno  step  and  four  leapfrog  steps,  with  a  time  in- 
crement of  six  minutes. 
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At  the  onset  the  pressure  levels  were  initialized  at 
standard  heights  and  temperatures  were  also  specified  by 
the  zonally  averaged  climatology.   The  model  was  started 
without  the  diabatic  forcing  function  until  a  quasi -steady 
circulation  was  established.   At  this  point  the  diabatic 
heating  was  introduced.   For  both  experiments  the  model 
runs  for  25  more  days  from  the  quasi-steady  state  with  the 
data  of  the  last  10  days  used  for  computing  the   time-mean 
fields.   All  data  presented  is  interpolated  to  constant 
pressure  surfaces  of  250  mb,  500  mb ,  and  850  mb. 
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III.   RESULTS 

The  principal  circulation  features  of  the  northern  summer 
monsoon  are  the  Tibetan  high,  the  Pacific  tropical  upper 
tropospheric  trough  (TUTT),  the  tropical  easterly  jet,  the 
Mexican  high  and  the  Atlantic  TUTT.   The  simulation  of  these 
features  by  the  model  are  compared  with  observations,  mainly 
those  by  Krishnamurti  and  Rodgers  (1970). 


A.   TIME  MEAN  FIELD 

1 .   Steady  Heating 

The  10  day  time-mean  fields  of  wind,  geopotential 
height,  temperature  and  vertical  velocity  are  shown  in  Fig- 
ures 4-7.   In  Figure  4a  the  250  mb  wind  field  shows  an  anti- 
cyclonic  wind  center  near  30°N,  52°E   which  resembles  the 
observed  Tibetan  high  although  the  position  of  the  generated 
feature  is  shifted  30°  to  the  west.   The  calculated  Pacific 
TUTT  extends  from  120°W, 40°N  to  180°W,  20°N,  approximately 
the  same  position  as  the  observed  TUTT.   The  TUTT  has  a 
northeast-southwest  tilt  and  the  amplitude  is  slightly  greater 
than  that  observed.   The  strongest  easterly  jet  is  developed 
by  the  model  in  the  southern  periphery  of  the  Tibetan  anti- 
cyclone.  The  jet  maximum  is  located  at  approximately  60°- 
90°E,  10°N,  with  a  speed  of  22  ms"   which  is  an  underestimate 
of  the  observed  maximum  speed  of  -30  ms"  . 

The  Mexican  anticyclone  is  centered  near  130°W,  25°N. 
This  feature  is  again  shifted  to  the  west  of  the  observed 
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position  by  25°.   The  amplitude  of  the  disturbance  compares 
well  with  the  observed. 

The  Atlantic  TUTT,  situated  in  the  mid-Atlantic  in 
the  area  of  15°N-40°N,  20°-60°W  is  well  developed  but  again 
the  amplitude  is  slightly  larger  than  that  observed.   As  for 
the  Pacific  TUTT,  there  is  no  apparent  phase  shift  problem 
when  compared  to  observations. 

There  are  no  good  observational  data  with  which  to 
compare  the  500  mb  flow  shown  in  Figure  4b.   The  850  mb  flow 
(Figure  4c)  clearly  shows  the  low-level  cyclonic  flow  below 
the  upper-level  Tibetan  and  Mexican  anticyclones.   Westward 
phase  shifts  of  circulation  centers  also  occur  at  low  level 
and  is  consistent  with  the  upper  level  flow. 

In  Figure  5a,  the  250  mb  geopotenti al  ,  the  Tibetan 
and  Mexican  highs  are  the  main  features.   The  Pacific  and 
Atlantic  TUTTs   are  weak  in  comparison  to  observations  but 
wel 1  def i  ned . 

Figure  6a,  the  250  mb  temperature  field,  shows  that 
the  Tibetan  and  Mexican  highs  are  warm  core  systems.   The 
TUTTs   in  both  the  Pacific  and  Atlantic  are  clearly  indi- 
cated by  cold  regions.   These  temperature  anomalies  are 
present  throughout  the  vertical  in  both  the  500  mb  field 
(Figure  6b)  and  the  850  mb  field  (Figure  6c). 

As  expected,  Figure  7,  which  shows  the  400  mb  verti- 
cal velocity  closely  resembles  the  horizontal  diabatic  heat- 
ing distribution.   This  is  consistent  with  the  warm  core 
energetics  of  the  monsoon  in  which  the  long  wave  kinetic 
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energy  is  converted  from  the  available  potential  energy 
generated  by  diabatic  heating. 

Figures  8  and  9  show  the  mean  divergence  and  vortic- 
ity  fields,  respectively.   The  250  mb  divergence  (Figure  8a) 
is  a  reasonable  reproduction  of  the  200  mb  divergence  field 
as  observed  by  Kri shnamurti .   This  indicates  that  the 
specified  heating  function  is  properly  chosen  for  simulating 
the  forcing.   The  vorticity  fields  are  representative  of 
the  wind  fields,  shown  in  Figure  4.   A  negative  vorticity  belt 
surrounds  the  globe  from  8°-25°N.   The  two  TUTTs  are  clearly 
shown  by  fingers  of  positive  vorticity  extending  into  the 
negative  vorticity  belts. 
2 .   Fluctuating  Heating 

The  10-day  mean  fields  of  the  fluctuating  heating 
forcing  function  are  Figures  10-15.   Only  the  250  mb  level 
of  all  the  fields  is  shown  except  in  the  cases  of  vertical 
velocity  (Figure  13)  and  the  time-mean  wind  fields  (Figures 
10a, b,c).   A  close  comparison  of  these  fields  with  the  time- 
mean  fields  for  the  steady  heating  case  could  not  reveal  any 
significant  differences  between  the  two  cases.   The  phase 
shifts  of  the  Tibetan  and  Mexican  highs  are  almost  exactly 
duplicated  by  the  fluctuating  heating,  resulting  in  almost 
identical  amplitudes  and  phases  as  those  for  the  steady 
heating. 

The  above  results  suggests  that  the  effect  of 
fluctuating  divergence  forcing  in  the  summer  monsoon  suggested 
by  Holton  and  Colton  is  not  prominent  enough  to  explain  the 
phase  problem  in  the  diagnostic  simulation  of  the  time-mean 
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field.   This  conclusion  is  true  if  the  amplitude  of  the  fluc- 
tuating component  in  our  model  is  representative  of  the 
actual  observed  amplitude.   Since  limited  observational  know- 
ledge (based  on  Kri shnamurti ,  1970)  does  indicate  that  the 
actual  amplitude  of  fluctuation  may  not  be  much  larger  than 
what  is  specified  here,  Holton  and  Colton's  hypothesis  seems 
unlikely  to  be  valid. 

B.   TRANSIENT  FIELD 

The  second  purpose  of  this  study  is  to  examine  the  synop- 
tic-scale transient  disturbances  associated  with  the  plane- 
tary-scale monsoon  circulation. 

1  .   Steady  Heating 

Figure  16  is  a  time-longitude  section  for  the  12-hour 
interval  200  mb  vorticity  field  at  12°N  from  day  15  to  day  25. 
In  the  western  hemisphere  from  0°-90°W  there  is  evidence  for 
strong  synoptic-scale  westward-propagating  wave  activity. 
The  wavelengths  of  these  disturbances  vary  from  3500  km  at 
30°W  to  2200  km  at  60°W.   In  some  cases  the  wavelength  can 
be  as  short  as  1000  km  at  180°W  to  90°E  where  the  disturbances 
appear  to  be  much  more  irregular  with  long  wavelengths  at  the 
beginning  of  the  time  series  and  sporatic,  short  wavelengths 
toward  the  end. 

Upstream  from  the  jet  maximum  near  50°E  the  disturb- 
ance activity  is  very  minimal  with  small  amplitude.   The  waves 
begin  to  grow  in  the  area  of  the  easterly  jet  maximum  and 
continue  to  grow  as  they  propagate  westward.   Downstream  of 
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the  easterly  jet  maximum  the  synoptic  waves  appear  to  be 
much  more  organized  comparing  to  upstream  of  the  jet. 

Krishnamurti  (1971)  and  Colton  (1972)  have  found 
that  the  upper-level  synoptic  waves  are  largely  due  to  the 
barotropic  instability  associated  with  the  easterly  jet. 
To  examine  this  possibility  as  the  generation  mechanism  of 
the  waves  observed  in  this  model,  we  plot  the  meridional 
gradient  of  the  basic  absolute  vorticity 

R    3  u 

where   3   is  the  gradient  of  the  coriolis  parameter  and   u 
is  the  time-mean  zonal  wind  component,  as  a  function  of  x 
and  y  in  Figure  17.   It  is  obvious  that  the  sign  of  this 
quantity  changes  within  the  north-south  channel  throughout 
the  domain.   Thus  it  is  quite  probable  that  the  barotropic 
instability  exists  everywhere  at  the  250-mb  level.   The  im- 
portant role  played  by  the  easterly  jet  in  this  regard  is 
reflected  by  the  fact  that  the  entire  latitude  band  of  12°N 
has  negative  values  for  this  quantity. 

The  result  that  the  wave  amplitude  is  minimum  up- 
stream of  the  easterly  jet  and  increases  to  a  maximum  down- 
stream further  suggests  that  the  mean  kinetic  energy  is  con- 
verted into  eddy  kinetic  energy  in  the  vicinity  of  the 
easterly  jet  maximum.   There  should  be  a  stable  area  upstream 
of  the  easterly  jet  maximum  where  the  amplitude  of  the  dis- 
turbances is  small.   The  amplitude  near  the  jet  maximum  may 
be  small  but  the  change  in  amplitude  is  large  as  the  wave 
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propagates  past  the  jet  indicating  that  the  jet  maximum  is 
the  most  unstable  region.   This  result  indicates  that  the 
most  active  area  for  synoptic-scale  wave  disturbances  is  in 
the  Atlantic,  Caribbean  and  Central  America,  not  in  the 
Pacific  and  Indian  Ocean  areas.   There  are  two  possible 
factors  in  the  model  that  influence  this:   (1)   The  model 
position  of  the  Tibetan  high  is  shifted  by  30°  to  the  west 
from  the  observed  position,  and  there  is  a  corresponding 
shift  in  the  jet  maximum  from  the  Indian  to  the  Arabian  Sea. 
The  most  active  wave  region  should  be  downstream  from  the 
jet  maximum  so  this  region  is  shifted  from  the  Arabian  Sea 
and  Africa  to  the  Atlantic  area.   Also  the  Atlantic  TUTT 
has  been  known  to  be  an  active  area  for  upper-level  cold  core 
lows  to  develop.   It  is  possible  that  the  synoptic  waves 
observed  here  are  also  a  reflection  of  the  cold  core  lows. 
On  the  other  hand,  the  absence  of  strong  disturbance  acti- 
vity in  the  Pacific  TUTT  region,  in  contrast  to  the  finding 
by  Colton  (1972),  may  also  be  explained  by  the  westward 
phase  of  the  easterly  jet  maximum.   (2)   This  model  contains 
a  horizontal  diffusion  which  tends  to  reduce  the  zonal  wind 
shear  and  in  so  doing  the  barotropic  instability.   This  dif- 
fusion was  not  used  by  Colton  (1972).   It  is  noteworthy, 
however,  that  even  in  this  model  the  barotropic  instability 
mechanism  is  still  found  to  be  prominent  downstream  from  the 
jet  maximum.   An  optimum  use  of  the  diffusion,  to  alleviate 
the  problem  associated  with  the  finite  differencing  and  at 
the  same  time  avoid  strong  artificial  smoothing  of  the  basic 
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wind  shear,  should  give  us  a  better  simulation  of  this 
mechanism. 

As  an  example  for  the  barotropi cal ly  unstable  wave, 
a  sequence  of  12-hour  interval  250  mb  wind  analysis  is  shown 
in  Figures  18  a-d,  from  0000  GMT,  day  22  to  1200  GMT  day  23. 
Here  we  shall  follow  the  development  of  three  waves  in  time. 
Their  positions  at  80°W,  4°-16°N;  60°W,  4°-16°N  and  30°W, 
4°-16°N  on  0000  GMT  day  22  are  indicated  by  dashed  lines. 
As  time  progresses  each  wave  propagates  to  the  west  with  a 
wavelength  of  -3500-4000  km.   The  meridional  tilt  of  the  de- 
picted waves  i  s  approximately  northwest-southeast.   From 
Figure  4a,  the  maximum  time-mean  easterly  zonal  wind  is  near 
12°N  in  the  vicinity  of  the  waves  in  question.   Thus  in  the 
12°-16°N  band  the  time-mean  zonal  wind  has  a  northeast-south- 
west shear  which  opposes  the  wave  tilt.   This  is  clearly 
another  indication  that  the  waves  are  barotropi cal ly  un- 
stabl e. 

2 .   Fluctuating  Heating 

Figure  20  shows  the  time-mean  meridional  gradient  of 
the  absolute  vorticity  for  the  fluctuating  heating  case  and, 
as  expected,  strongly  resembles  Figure  17  for  the  steady  heat- 
ing case.   There  is  no  question  on  the  plausibility  for  baro- 
tropic  instability  to  exist. 

The  time-section  of  the  mean  vorticity,  Figure  19, 
also  points  out  that  the  greatest  wave  activity  is  downwind 
of  the  maximum  jet.   The  amplitude  of  the  waves,  however, 
appears  to  be  greater  in  this  fluctuating  heating  case,  though 
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the  range  of  wavelengths  remains  the  same  as  the  steady  heat- 
ing case.   There  is  also  no  noticeable  change  in  phase  speed 
of  ~7.6°  day"  .   This  increase  in  amplitude  is  readily 
apparent  in  the  waves  generated  in  the  Atlantic  area  of  20°- 
90°W,  8°-16°N.   This  result  indicates  that  the  fluctuation 
of  the  total  forcing  field  may  actually  enhance  the  barotropic 
energy  conversion.   In  Figure  21,  by  following  the  waves 
initially  located  at  85°W,  8°-16°N;  55°W,  10°-16°N  and  30°W, 
8°-16°N;  it  is  found  that  the  tilt  of  the  waves,  northwest- 
southeast,  again  opposes  the  horizontal  shear  of  the  basic 
wind . 
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IV.   CONCLUSIONS 


One  major  conclusion  of  this  study  is  that  the  westward 
phase  shift  problem  of  the  simulated  Tibetan  and  Mexican 
highs  in  the  absence  of  very  strong  damping  is  not  likely  to 
be  explainable  by  the  influence  of  the  fluctuating  heating. 
On  the  other  hand,  the  tropical  upper  tropospheric  troughs 
do  not  exhibit  this  phase  problem.   This  may  indicate  the 
importance  of  damping  due  to  cumulus  transport  of  momentum 
and  vorticity  because  cumulus  convection  is  predominant  below 
upper  level  high  pressure  areas  only. 

In  this  model  the  upper-level  barotropic  instability 
seems  to  be  quite  important  for  most  of  the  area,  even  though 
the  easterly  jet  was  not  developed  to  the  fully  observed 
strength,  perhaps  due  to  the  inclusion  of  the  somewhat 
strong  horizontal  diffusion.   This  instability  seems  par- 
ticularly important  when  time-fluctuations  of  the  planetary- 
scale  is  included  through  the  imposed  fluctuating  heating. 
Thus  the  local  barotropic  instability  at  upper  levels  may 
be  more  important  than  that  implied  by  the  horizontal  shear 
of  the  time-mean  wind.   Further  investigations  on  this 
mechanism  would  be  desirable  to  see  whether  the  observed 
upper  level  disturbances,  such  as  the  cold  core  lows,  are  of 
this  nature.   Incorporation  of  cumulus  parameterization  and 
topography  in  the  model  should  further  enable  us  to  study 
the  possible  implication  of  this  upper  level  barotropic  in- 
stability in  the  vertical  development  of  other  tropical 

disturbances. 
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In  addition  to  the  phase  shift  problem,  Holton  and 
Colton  (1972)  also  found  yery    large  amplitudes  for  long 
waves  in  the  case  of  weak  damping  in  their  linear  model. 
Our  results  do  not  have  this  amplitude  discrepancy.   This 
indicates  that  the  transient  waves  due  to  non-linear  baro- 
tropic  energy  conversion  serves  as  an  energy  sink  of  the 
planetary-scale  flow. 
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Figure    2.      Vertical    levels:       Sigma    is    the   vertical    coordinate. 
V   is   the   horizontal   vector   velocity.      tt    is    the    terrain   pressure; 
T    is    the   temperature   and   $    is    the   geopotential . 
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Figure  16.   250  mb  vorticity  time-longitude  section  at 

12°N,  steady  heating.   Hatched  area  indicates 
positive  values. 
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Figure  19. 


250  mb  vorticity  time-longitude  section  at 
1  2  °  N  ,  fluctuating  heating.   Hatched  area 
indicates  positive  values. 
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